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ABSTRACT 

By selecting astrometric and photometric data from the Sloan Digital Sky Survey (SDSS), the Lepine 
& Shara Proper Motion North Catalog (LSPM- North), the Two Micron All Sky Survey (2MASS), 
and the USNO-Bl.O catalog, we use a succession of methods to isolate white dwarf candidates for 
follow-up spectroscopy. Our methods include: reduced proper motion diagram cuts, color cuts, and 
atmospheric model adherence. We present spectroscopy of 26 white dwarfs obtained from the CTIO 
4m and APO 3.5m telescopes. Additionally, we confirm 28 white dwarfs with spectra available in the 
SDSS DR7 database but unpublished elsewhere, presenting a total of 54 WDs. We label one of these 
as a recovered WD while the remaining 53 are new discoveries. We determine physical parameters 
and estimate distances based on atmospheric model analyses. Three new white dwarfs are modeled to 
lie within 25 pc. Two additional white dwarfs are confirmed to be metal-polluted (DAZ). Follow-up 
time series photometry confirms another object to be a pulsating ZZ Ccti white dwarf. 
Subject headings: Catalogs - Proper motions - Stars: variables: general - Surveys - white dwarfs 



1. INTRODUCTION 

White dwarfs (WDs) are of interest in a variety of sub- 
fields in astrophysics because of their unique ability to 
act as cosmic chronometers. Theoretical cooling mod- 
els provide a means for dating WDs from photometry 
alone, and this process has provided constraints on the 
age of the Ga lactic disk using cool WD samples (e.g., 
iLeggett et al.l fl998) . Because the coolest (and hence old- 
est) WDs are less luminous, the nearby representatives 
provide the best opportunities for accurate characteriza- 
tion. 

The present under standing of the local WD sample is 
somewhat uncertain. iHolberg et al.l (|2008f ) estimate that 
the WD sample is ~80% complete to 20 pc. Includ- 
ing only WDs w i th ac curate trigonometric parallaxes, 
iSubasavage et al.l (|2009l ) conclude that the sample is only 
~47% complete to 25 pc. The vast majority of the in- 
completeness arises from the coolest WDs that remain to 
be discovered. Additional discoveries of nearby WDs will 
strengthen the completeness statistics and provide valu- 
able model parameter constraints through trigonometric 
parallax measurements. 

The Sloan Digital Sky Survey f SDSS: r^AeFan[2000l) 
has been exceptional at identifying WDs. The lastest re- 
lease as of the writing of t his manuscript. Data Release 
8 (DR8: lAihara et al.ll2011| ). reports sky coverage of over 
14,500 deg^. Imaging data are collected in the ugriz 
bands with a 50% completeness limit at r = 22.5 for point 
sources, and spectra have been observed for over half of a 



|csayres@u. Washington. edul 

^ Department of Astronomy, University of Washington, Box 
351580, Seattle, WA 98195 

^ US Naval Observatory, 10391 West Naval Observatory Road, 
Flagstaff, AZ 86001-8521 

^ Departement de Physique, Universite de Montreal, C.P.6128, 
Succursale Centre- Ville, Montreal, QC H3C 3J7, Canada 

* Visiting Astronomer, Cerro Tololo Inter-American Observa- 
tory. CTIO is operated by AURA, Inc. under contract to the 
National Science Foundation. 

5 CTIO REU 2010. 



million stars. Prior data releases from SDSS have led to a 
proliferation of WD discoveries. lEisenstein et al.l (|2006( ) 
(hereafter SDSS-E06) roughly doubled the number of 
previous, spectroscopically id entified WDs using spec- 
tra fr om the DR4 database (jAdelman-McCarthv et aLl 
|2006[ ). though most were hotter than ~7000 K because 
of biases in the SDSS spectrosco pic target s electio n pro- 
cess. The observational efforts of lKilic et all (|2006l . l2010[) 
led to spectroscopic confirmation of more than 100 cool 
WDs, where target selection was based on photometry 
and proper motions con tained in the SDSS DR7 database 
(jAbazaiian et al.ll2009l ). 

In this work, we utilize data from a suite of sur- 
veys and catalogs, in cluding SDSS, the Tw o Micron All 
Sky Survey (2MASS : ISkrutskie et "all[2006l) . USNO-Bl.O 
(iMonet et al.l l2003f). a nd the proper motion survey of 
iLepine fc Sharal ()2005| ) to identify new WD candidates. 
From these data, we prioritized targets for follow-up 
spectroscopy, with emphasis on cool, nearby WDs as 
well as hotter, and potentially variable (i.e., ZZ Ceti) 
WDs. We present spectra and physical parameters, de- 
rived from spectral energy distribution (SED) and model 
atmosphere analyses, for 25 newly discovered WDs and 
one recovered WD (see 14. 5|) . In addition, we present 
physical parameters for 28 new WDs that have been 
spectroscopically confirmed by SDSS DR7 spectra but 
are unpublished elsewhere. 

2. TARGET SELECTION 

2.1. Methodology 

Target selection was based on a combination of con- 
straints applied to optical and near-IR photometry, 
proper motion, and mo del- adherence. Two independent 
target selection passes were employed due to high con- 
tamination rates experienced with the initial approach. 
We first queried the SDSS DR7 database for objects 
with the following stipulations: proper motion > 100 
mas yr~^, g magnitude < 19.5, and declination < -f30°. 
The declination constraint was used to ensure all tar- 
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gets could be observed from Cerro Tololo Inter- American 
Observatory (CTIO), where the fi rst round o f spec - 
troscopic observations was taken. iMunn et al.l (|2004[ ) 
have joined the USNO-Bl.O Catalog with SDSS astrom- 
etry (available in the propermotions table of the SDSS 
database), which allowed us to query for high proper 
motion objects directly from SDSS. In a second pass, 
w e began with proper rn otion objects from the catalog 
of lLepine fc Sharal (J2005D containing objects with proper 
motions >150 mas yr""'^. For these objects, we required 
a match in SDSS and employed the g magnitude < 19.5 
limit with no declination constraint. In both approaches, 
candidate objects were required to have near-IR photom- 
etry from 2MASS. Cross- matching between catalogs was 
achieved using the table matching routines of TOPCATO 
([Taylor 2005). 

The technique of using reduced proper motion ( RPM) 
to id entify WDs has been used for decades (e.g., iJonesI 
Il972| ). Conceptually, RPM is used as a proxy for ab- 
solute magnitude. The two quantities are connected by 
the inclusion of tangential velocity. Thus, proper mo- 
tion coupled with color and apparent magnitude serve 
to separate generally blue, lower luminosity, and larger 
tangential velocity WDs from low-metallicity halo subd- 
warfs (SDs) and main-sequence stars. 

Reduced proper motion is defined as 

i7„ = m + 5 log ^ + 5 = M + 5 log Vtan - 3.379 (1) 

where m is apparent magnitude, M is absolute magni- 
tude, n is proper motion in arcseconds yr~^, and Vtan 
is tangental velocity in km s^^. In this study, only tar- 
gets satisfying Hg > 15.136 -f 2.727{g — i) were kept as 
WD candi dates, a reduced proper motion diagram cut 
defined bv lKilic et al.l ([2006) that eliminates most of the 
low-metallicity SDs and virtually all main-sequence stars. 
A photometric color cut J — Ks < 0.5 was employed 
to better constrain initially large sam ples of WD candi- 
dates. This color cut was adopted from lSubasavage et al.l 
([2008) where they found this to be a clear delimiter be- 
tween WDs and SD contaminants. All candidates were 
verified to have noticeable proper motions. These by-eye 
verifications were carried out by blinking between digi- 
tized POSS I and POSS II epoch s for each target usin g 
the Aladin interactive sky atlaqH (jBonnarel et al.l l2000V 
One final model-adherence step was implemented to 
better probe the cooler regime (Toff < 7000 K) that is 
plagued by subdwarf contaminants. Remaining targets 
were crudely compared with WD synthetic photometry 
models to estimate Toff values and distances prior to ob- 
servations. Targets poorly represented by the models as 
defined below were discarded. We used pure-hydrogen 
atmospheric WD modelqj and assumed log g = 8.0. The 
model grid was cubic spline interpolated to give a tem- 
perature resolution of 10 K. Both modeled and candidate 
photometry were converted to F\ flux v alues by the pre- 
scription of !Holb erg fc Bergeronl ()2006[ ). The resulting 
SEDs were normalized by the r band flux value. A best 
fit was determined using a chi-square minimization be- 
tween model and target normalized flux values. If the 
photometric error was greater than 0.1 mag then that 
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passband was ignored for fitting purposes - this only oc- 
curred in the near-IR passbands as the SDSS magnitudes 
are significantly better down to our adopted magnitude 
limit oi g = 19.5. A p value for goodness-of-fit was de- 
termined using common chi-square lookup tables. We 
elected to keep only WD candidates with p > 0.95. 

The fitting process provided estimates of T^s and dis- 
tance for each target in the final sample. After spectro- 
scopic observations, all bona fide WDs were modeled in 
a more robust fashion for physical parameter determina- 
tions (described in 14.11 and presented here) and should 
not be confused with this pre-observation model fitting. 
Finally, all previously identified objects were discarded 
from our target list, the vast majority of these were pub- 
hshed WDs. 

2.2. Completeness Estimates 

We estimate our target-selection completeness by look- 
ing at the recovered fraction of SDSS-E06 WDs that met 
our criteria for both samples (Pass 1 and 2). For a homo- 
geneous comparison, all of the relevant DR8 data (i.e., 
RA, Dec, proper motion, ugriz, JHKs) were extracted 
for the SDSS-E06 sample by matching platc/fiber/mjd 
designations. A total of 39 objects out of 9316 in SDSS- 
E06 were not recovered because the spectra are not in- 
cluded in the DR8 database (hence, no plate/fiber/mjd 
designation). We confirm that they are available in the 
DR7 database. Nevertheless, these few missing objects 
should not significantly affect the completion statistics. 

To compare the SDSS-E06 sample (with DR8 data) 
to our Pass 1 sample, we implemented identical obser- 
vational criteria (i.e., proper motion > 100 mas yr~^, 
g < 19.5, Dccl. < +30°, J - Kg < 0.5) to the SDSS- 
E06 sample. Our Pass 1 sample was then positionally 
cross-matched to this sample and 19 out of 21 were re- 
covered, implying a completeness of 90%. To estimate 
the completeness of our Pass 2 sample, the SDSS-E06 
sample (with DR8 data) was positionally cross-matched 
to the LSPM catalog. For all common objects, identi- 
cal observational criteria (i.e., g < 19.5, J — Ks < 0.5) 
as implemented in our Pass 2 were applied. Our Pass 
2 sample was compared via LSPM name to the remain- 
ing objects once these criteria were applied and 31 out 
of 35 were recovered, implying a completeness of 89%. 
Note that in both cases, our adopted limiting magnitude 
excluded the vast majority of WDs in the SDSS-E06 cat- 
alog thereby leaving us with fairly small numbers with 
which to estimate completeness. 

To better understand how our selection criteria affects 
our completeness, the cuts relating to direct observables 
(i.e., proper motion, magnitude, color, and declination) 
were performed first. These samples represent the 100% 
complete subsample of SDSS-E06 relevant for each of our 
passes. We then implement the RPM cut and, in both 
cases, recover 100% of the subsample. It was only in the 
model-adherence cut that recovery rates fell below 100%. 
Because only a total of 5 unique WDs were not recovered 
(one unrecovered object belonged to both samples), we 
looked more carefully at each one to ensure our model- 
adherence criteria were robust. WD 0756-1- 437 is a mag- 
netic DA with a field strength of --300 MG ([Kiilebi et al.l 
I2009f ). WD 1156-^132 is classified as DQpec with deep 
swan-band carbon absorption. WD131 1-1-129 is classified 
as a DBA and has the largest hydrogen abundance of the 
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iBergeron et~all (|2011[ ) sample of DBAs [log (H/He) = - 
2.90]. The remaining two WDs, WD 1559+534 and SDSS 
J101436.01+422622.0 appear to be normal DA white 
dwarfs that only marginally failed our model-adherence 
criteria with p > 0.9. It stands to reason that the first 
three would fail our mo del- adherence criteria because 
their SEDs are not well-represented by a pure-H model 
atmosphere. Thus, while our completeness suffered be- 
cause of our model-adherence criteria, the survey was 
better able to successfully probe the parameter spaces 
dominated by contaminants (Tcff < 7000 K) to identify 
nearby WDs as discussed in 14.41 



3. DATA AND OBSERVATIONS 
3.1. Astrometry and Nomenclature 

For new discoveries, we determine WD names 
in the conventional manner (fully described in 
iSubasavage et al.ll2007() by using the target's epoch 1950 
equinox 1950 coordinates. 

As stated previously, proper motions were initially 
extracted from the SDSS DR7 database and con- 
sisted of combined U SNO-Bl.O pl u s SDS S astrometry 
(|Munn et all I2004D . iMunn et al.l (|2004D demonstrate 
that the proper motions derived from combined astrome- 
try improves those contained within the USNO-Bl.O cat- 
alog by ^25% when compared to bright, n on-moving, 
spectroscopically-confirmed QSOs. However, iKilic et aD 
(pbOQl found that high proper motion objects with neigh- 
boring sources within 7" were more likely to have incor- 
rectly measured p roper motio n s. To remove this source 
of contamination IKilic et ahl ()2006[ ) discarded objects 
with neighbors within 7". We did not implement this 
criterion to avoid the possibility of missing true WDs. 

After an initial night of spectroscopic observations at 
CTIO (discussed in Section [?75)) . we realized 14 SD and 
main-sequence contaminates with incorrect proper mo- 
tion values from the SDSS query. By-eye verification con- 
firmed these to be bona fide proper motion objects, but 
proper motion magnitudes and position angles were er- 
roneous, suggesting mismatches when combining USNO- 
Bl.O and SDSS astrometry. To mitigate this effect, we 
initi ated a seco nd query using the LSPM-North Cata- 
log (jLepine fc S hara 2005) as a starting point for proper 
motion values and cross-matched with the SDSS DR7 
and 2MASS databases. LSPM-North objects were not 
initially verified by eye for proper motion confirmation 
as this was painstakingly done by the original authors. 
This lead to a significant reduction in contaminants. 

Following observations, proper motions for all spec- 
troscopic confirmations (WDs and contaminants) were 
doubl e checked using the SuperCOSMOS Sky Survey 
(SSS; IHamblvet an[200l . In cases where the USNO- 
Bl.O-l-SDSS and SSS proper motions were discrepant, 
a by-eye inspection was performed to confirm position 
angle. For these objects, we adopt the proper motions 
extracted from the SSS. 

WD names, alternate names (LSPM-North where 
available and SDSS otherwise), epoch 2000 coordinates 
and adopted proper motions are listed in Table 1 for 
both WDs {top) and contaminants (less the WD name, 
bottom). 

3.2. Photometry 



The psfMag values in the optical ugriz passbands were 
extracted from the SDSS DR8. The target selections 
were performed using DR7 photometry but in the in- 
terim, DR8 was released and thus, we use these values 
for SED modeling discussed in Section |4?T] We list these 
values and their corresponding errors in Table 1. 

As a requirement, all candidate WDs had to contain 
JHKs data in 2MASS. Given that the majority of these 
targets are near the faint limit of 2MASS, we utilized 
the UKIRT Infrared Deep Sky Survey (UKIDSS) DR6 
Large Area Survey to supplement JHKs for new discov- 
eries where avai l able. The UKIDSS project is defined in 
ILawrenceetaH ([Ml . U KIDSS uses the U KIRT Wide 
Field Camera (WFCAM; ICasah et al.l 120071) The pho- 
tometric system is described in'Hewett et al.' ('2006'), and 
the calibration is described in Hodgkin ct al. (2009). The 
pipeline processing and science archive are descri bed in 
Irwin et al. (2009, in preparation) and iHamblv et al.l 
(2008). UKIDSS magnitudes were tran sformed to the 
2MASS system using the prescription of iHodgkin et al.l 
(2009;). The JHKs values and corresponding errors are 
listed in Table 1. Converted UKIDSS photometry is 
listed whenever errors in JHKs are less than 0.05. Oth- 
erwise, the values arc from 2MASS. 

3.3. Spectroscopy 

Prior to observations, 28 previously unidentified ob- 
jects had spectra in the SDSS DR7 database confirming 
their WD nature. These objects are noted in Table 1, 
but we do not include their spectra as they are freely 
available in the SDSS database. 

Spectroscopic observations were conducted throughout 
2010 and early 2011 from both CTIO and Apache Point 
Observatory (APO). At CTIO the 4m Blanco Telescope 
and Ritchey-Chretien Spectrograph were used. We se- 
lected the KPGL3 grating that covered a wavelength 
range of 3600-7000 A. Observations in 2010 were taken 
using a 2'.'0 slit width, oriented due north-south, to pro- 
vide spectral resolution of 6 A. Observations in 2011 
were taken using a 4'.'0 slit width to minimize light loss 
from differential refraction as the slit was not rotated 
to the parallactic angle but rather was kept fixed again 
due north-south. With this configuration, the resolu- 
tion degraded slightly to 8 A. As can be seen in Figure 
[Jl which contain spectra from both runs, differential re- 
fraction was not a significant problem during the first 
run, nor was the loss of resolution during the second run 
detrimental for the purpose of spectral classification. 

At APO, the Dual Imaging Spectrograph (DIS) on the 
ARC 3.5m telescope was used with the B400+R300 grat- 
ings for spectroscopic follow-up. We achieved 6 A reso- 
lution spanning 3350-9260 A, but fringing became prob- 
lematic beyond '^7000 A. Red and Blue spectroscopic 
channels were reduced independently. Red channel flux 
near the dichroic at approximately 5400 A was found 
to be highly variable, so wavelengths less than 6300 A 
were omitted in the red spectral channel. The default 
slit width of 1'.'5 was used for the DIS spectrograph. 

During all of the spectroscopic observing runs, flux 
standards were observed each night for flux calibration 
and HeNeAr lamps were taken at each telescope pointing 
for wavelength calibration. Two spectra of each target 
were obtained to permit cosmic ray rejection. Data were 
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reduced using standard IRAlfl routines. 

We confirm 19 new WDs from CTIO 4m observations. 
Figure [1] presents 14 DA discoveries and one DA recov- 
ery (see Section |475|) . Figure [2] displays one DQ and two 
DC WDs. Figure [3] shows spectra and model fits for two 
DAZ WDs with detectable Ca II absorption features at 
3933 and 3968 A. From CTIO, 18 contaminants were ob- 
served and their spectra are plotted in Figure|31 Defining 
absorption features in these contaminant spectra are due 
to metal and molecular content (Ca II, CH, MgH). 

We confirm 6 new WDs from APO 3.5m observations: 
three DA, two DC and one DZ with Ca II absorption 
in an otherwise featureless spectrum. The upper panels 
of Figure [SI show these discoveries. We experienced only 
two contaminants (both likely SDs) and their spectra are 
plotted in the lower two panels of Figure [5] 

4. ANALYSIS 

4.1. Modeling of Physical Parameters 

Our model atmosphe res for WDs a r e simi lar to those 
described at lengt h in iLiebert et al.l (I2005L and refer- 
ences therein) and iBergeron et aL (ll995aD. with s evera l 
imp rovements discussed iii iTremblav fc BergeronI (|2007[ ) 
and iTremblav et "all ()2011[ ). In particular, we now make 
use of the improved Stark broadening profiles for the hy- 
drogen lines developed by Tremblay & Bergeron (2009). 
Our models for DQ and DZ stars, which include metals 
and molecules in the equation of st ate and opac i ty cal - 
cula tions, are des c ribed in detail in lDufour et al.l ()2005l ) 
and lDufour et al.l ()2007[ ). respectively. 

Table 2 contains SED-derived Toff (column 2) and 
distance (column 4) for WDs with pure H or pure 
He atmospheres (denoted in column 3). A complete 
di scussion of our SED-fi tting procedure can be found 
in IBergeron et "all (pOOl . For "polluted" WDs (DAZ, 
DQ, DZ), these parameters rely on the iterative, com- 
bined SED a nd spectral fitting procedures defined in 
iDufour et all (2005, 2001). Metal abundances derived 
from the model fits for these targets are listed in the 
notes section of Table 2. For all cases, we assume log g = 
8.0. Spectral subtypes (column 5) are d etermined for the 
DA WDs using the temperature index of lMcCook fc SionI 
([T999f), where the temperature index equals 50,400/TefT 
. In addition, spectro scopic line profi l e fitti ng was per- 
formed as described in iBergeron et al.l (|1992| ) for all DAs 
(and one DBA - WD 1457-^249) with sufficient line ab- 
sorption and signal-to-noise to produce a reliable fit. 
These results are listed in columns 6 and 7 of Table 2. 

4.2. Metal-Polluted DA White Dwarfs 

Figure [3] displays spectra and model fits for WD 
0920-^012 and WD 1408-^029. These two DA WDs ex- 
hibit Ca II H & K features and earn the classification 
of DAZ. Hypotheses to explain these spectral features 
include enrichment from (1) the interstellar medium or, 
(2) debr is disk accretion, with the latter being heavily 
favored (jFarihi et al.llMol ). 

WD 0920-^012 is estimated to be at 33.4 pc, while 
WD 1408-1-029 is estimated to be at 26.5 pc, and both 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatry, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 



are modeled to have a log (Ca)/(H) — -9.0. Given their 
relatively bright apparent magnitudes (g ^ 17), these 
targets would be excellent candidates for follow-up IR 
studies to possibly detect emission from the accretion 
disk and better characterize the system. 

4.3. ZZ Ceti White Dwarfs 

Two of our confirmed WDs were modeled to lie in 
the ZZ Ceti instability strip: WD 1419-1-062 and WD 
2102-1-233. Figure E| shows the instability strip as re- 
cently redefined by IGianninas et"aLl (|2011| ) using im- 
proved model spectra with the new Stark profiles de- 
scribed above. 

Differential photometry was performed on both ZZ 
Ceti candidates. Candidate flux was normalized by 
bright non-varying stars in the same field to obtain a 
differential light curve. The light curve was then nor- 
malized by the mean value to obtain fractional variation 
about the mean. Frequency content was analyzed using 
the magnitude of the Fourier Transform (FT) of the time 
series data. 

WD 2102-^233 was observed using the CTIO 0.9m tele- 
scope using the fuh 13f6 field and the BG 40 fiher. The 
time cadence was ~50 seconds. As can be seen in Fig- 
ure [71 we identify a dominant pulsational period at ^^800 
seconds with an amplitude of ^^2.6%. 

WD 1419-1-062 was observed using the CTIO 1.0m tele- 
scope, with a 20.'0 field, and the BG 40 filter. We detect 
no obvious pulsations from these observations (see Fig- 
urejS]). From the FT, the observed noise level is ^0.25%, 
with no obvious peaks. With our less-than-optimal sam- 
pling rate (^100 seconds), we could only hope to detect 
periods > 200 seconds based on the Nyquist sampling 
theorem, so this object may be pulsating outside our de- 
tection sensitivity. However, considering its location in 
Figure [Bj lying near the red edge of the instability strip, 
we would expect a longer period variable. 

4.4. Nearby White Dwarfs 

We find three new WDs with distance estimates within 
the 25 pc horizon of interest (WD 1338-1-052 at 13.7 ± 2.7 
pc, WD 1630-H089 at 13.2 ± 2.3 pc, and WD 2119-^040 
at 22.1 ± 3.6 pc) that has been adopted from the Cat- 
alog of Nearby Stars (CNS; [G licsc & Jahreifi 199 1]) and 
the NStars Database (|Henrv et al.ll2003l ). iHolberg et al.l 
(|2008) determine that the 20 pc WD sample is -80% 
complete based on the assumption that the 13 pc WD 
sample is complete. Two of these WDs, and possibly the 
third if it is slightly less distant than expected, will lie 
within 20 pc if proximity is confirmed. Moreover, WD 
1338-f 052 and WD 1630+089 are modeled to lie at a dis- 
tance very near to 13 pc. If either object proves to lie 
within 13 pc, the local WD population must be denser 
than previously thought. With any new WD discover- 
ies within 13 pc, a constant-density extrapolation would 
increase the amount of missing WDs in the 13 pc to 20 
pc range. The exact current completeness statistics will 
depend on more robust distance determinations. 

These three objects are being observed for trigonomet- 
ric parallaxes via the Cerro Tololo Inter- Americ an Obser- 
vatory Parallax Investigation (CTIOPI: pao et al. 20051 
'Henrv et al." '2006; 'Sub asavage et al.l 120091: iRiedel et atl 
,2010; iJao et al.,i2011.) program to confirm proximity. To 
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the best of our knowledge, the two WDs estiniated to be 
within 20pc are the nearest WDs discovered using SDSS 
data, if trigonometric parallaxes confirm proximity. Con- 
tinuing discoveries of the coolest WDs, especially those 
nearby and thus suitable for parallax measurements, will 
provide anchor points for WD atmospheric models that 
are vital to widely-used cosmic dating techniques. 

4.5. Notes in Individual Systems 

WD 0351-002 is also known as SA 95- 42 a. n d is a 
popular spcctrophotometric standard star. lOkd ()199Clf ) 
lists this object along with 24 other spcctrophotomet- 
ric standard candidates, including spectra. For reasons 
unclear, this is one of thr ee ob j ects t hat do not have a 
spectral type in Table 1 ofl Qkd (|1990( ) yet the spectrum 
published in that same work shows broad Balmer absorp- 
tion indicative of a DA WD. We include it here as a new 
discovery. 

WD 0412+065 is also kn own as GD 59. It was 

classified as a WD suspect by iGiclas et al.l ()1965( ) , but 
we found no spectroscopic confirmation in the literature 
even though it is listed as a WD in Simbad. We include 
it here as a new discovery. 

WD 1402+065 is also known as PC 1402+065 and 
was previou sly labeled as a subdwarf in the catalogs of 
IGreen et al.l ((1986) and subsequently i ncluded in the sub- 
dwarf catalog of lKilkennv et al.l ()1988l ) , based on the pre- 
vious determination. Our spectrum (see Figure[T]) clearly 
shows broad Balmer absorption and thus, we include it 
here as a new WD discovery. In fact, this WD is the 
hottest of those spectroscopically observed from CTIO 
plotted in Figure [I](Teff = 26,190 K). 

WD 1419+062 is also known a s PG 1419+062 and 
was first pu blished by [G reen (1980) as a DA. However, 
it appears in IGreen et al. (1986) with a "sd" designation 
indicative of being a hot subdwa rf. It then appears in the 
catalog of lKilkennv et al.1 (I1988D. aga in classified as a hot 
subdwarf based on the iGreen et aD (1986) designation. 
Here we confirm it to be a bona fide DA WD and suspect 
thelGrecn ct al. (1986) designation is a typo. We include 
the spectroscopy and the modeled parameters yet do not 
classify this object as a new discovery. This object is also 
discussed in Section l473l as a ZZ Geti candidate. 

WD 1434+159 is also kn own as GD 168. I t was 
classified as a WD suspect bv IGiclas et al.l ()1965| ). but 
we found no spectroscopic confirmation in the literature 
even though it is listed as a WD in Simbad. We include 
it here as a new discovery. 

WD 1457+249 has a spectrum dominated by He and 
was SED-modeled using a pure He model atmosphere. 
During the spectroscopic line fitting analysis, it became 
evident that, because of weak Balmer H/3 absorption, 
trace amounts of H existed in the atmosphere. Thus, 
this object is classified as a DBA with the best-fitting 
model including log (H/He) = -5.8 ± 0.1. 

5. DISCUSSION 

We present and characterize 54 WDs. Of these, 
26 objects were spectroscopically observed from either 
CTIO or APO (including recovered known object WD 
1419+062, as discussed in Section HTS)) . while the remain- 
ing 28 objects have spectra available in the SDSS DR7 
database for a total of 53 new WDs. It is likely the SDSS 
DR7 spectra will be contained in the forthcoming DR7 



White Dwarf Catalog described by iKleinmanI (|2010D . 
Additionally, 21 objects overlap (17 from SDSS spectra 
and 4 photometrically-s elected) between thi s study and 
a recent publication by iGirven et al.l (|2011h to identify 
DA WDs in SDSS. 

We find three WDs expected to be in the local neigh- 
borhood, with the closest (WD1630+089) expected to be 
13.2 pc distant. Twenty-one WDs are modeled to have 
Tcff < 7000 K, with four of those objects modeled to have 
Teff < 5000 K. We also confirm WD 2102+233 as a new 
ZZ Ceti star. 

Our survey is summarized in Table 3. Briefiy, our se- 
lection criteria were designed to optimally identify cooler 
WDs in parameter spaces where contaminants dominate. 
As such, we realized some contamination in our final 
sample, both identified by our spectra and also from 
the literature (Table 3 column 10) as listed by Simbad. 
Our methodology was initially impacted by problematic 
proper motions leading to significant contamination by 
subdwarfs and main-sequence stars (Pass 1 in Table 3). 
Once correct proper motions are applied to these objects 
(see Figure JH]) , it is clear the contamination would have 
been largely avoided as these objects fall above our RPM 
cut. Our corrected c ontamination is generally consistent 
with iKilic et all (|2010. ) who show decreasing contamina- 
tion rates with increasing Vtan cutoffs. They quote a 1.3% 
contamination rate for T4an ^ 30 km s~^ and indeed, we 
find one contaminant scattered within the Vtan ^ 30 km 
s^^ cutoff. While we show yet another example of the 
effectiveness of RPM to isolate WDs, there is an inherent 
bias against the slowest moving WDs. It is unlikely that 
this bias will be completely removed until a magnitude- 
limited astrometric survey is conducted, such as Gaia or 
LSST. 

Only a subset of our WD candidates were able to 
be observed because of telescope time constraints. The 
last column in Table 3 shows the number of targets left 
unobserved. Furthermore, all of our targets were con- 
strained to the SDSS footprint thereby covering only a 
fraction of the sky (Table 3 column 3) that is largely 
weighted towards northern declinations. For these rea- 
sons, we expect that more WDs of interest remain undis- 
covered. Highly anticipated photometric surveys such as 
Pan-STARRS and the Dark Energy Survey (DES) will 
certainly aid in the discovery of additional WDs. With 
this in mind, we publish our observed contaminates to 
add to those already known to have similar p hotometric 
properties a nd proper motions to WDs (e.g.. IKilic et al.l 
120061 I2010D . By identifying large samples of contami- 
nants now, disentangling them from WD candidates may 
be easier in the future u sing e mpirical and statistical 
methods. :Newling et al] ()2011h provide an example of 
statistical techniques based on photometry to classify su- 
pernovae. For similar techniques to work in our context, 
large robust training sets of both SD and WD exemplars 
are needed before any statistical assertions can be made, 
and applied to large data sets. 
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SDSS 


Null 


CTIO 


0.027 


SDSS 


Null 


CTIO 


Null 


SDSS 



CO 






d 



Table 1 — Continued 



Nar 



Nar 



1516+128 
1519+141 
1531+024 
1547+481 
1618+208 
1630+089 
2051-051 
2102+233 
2119+040 
2144+156 



SDSS J151911.99+123946.5 
LSPM J1521+1358 
SDSS J153417. 49+021848.0 
LSPM J1549+4a02 
LSPM J1620+2044 
LSPM J1632+0851 
SDSS J205342. 94-045939. 8 
SDSS J2 10452. 71 +233320. 5 
LSPM J2122+0413 
LSPM J2146+1550 



15 19 11. 
15 21 20. 
15 34 17. 

15 49 27. 

16 20 27. 
16 32 33. 

20 53 42. 

21 04 52. 
21 22 12. 
21 46 32. 



98 +13 58 

50 +02 18 

16 +48 02 

73 +20 44 
18 +08 51 
95 -04 59 

74 +23 33 

35 +04 13 

36 +15 50 




3 


SDS 


,t 


LEP 





SDS 


3 


LEP 


R 


LEP 




SDS 


fi 


SDS 




SDS 


4 


SDS 


4 


LEP 



17 


613 





012 


7 


169 






t 


008 


6 


680 


16 


680 


t 


007 


6 


318 


IS 


053 


C 


015 


7 


403 




233 


t 




6 




16 


445 


t 


007 


5 


345 






t 




5 


512 


16 


238 


t 


006 


5 


889 






t 






200 


17 


137 





010 


6 


603 



198 


005 17 


277 





006 17 


459 


788 C 


005 16 


937 




006 17 


116 


272 C 


004 16 


302 


L 


005 16 


416 


154 C 


005 17 


087 


C 


006 17 


121 


734 C 


004 16 






005 16 


857 


877 C 


005 14 


709 


L 


005 14 


353 


632 C 


004 15 






004 15 


942 


036 C 


003 16 


148 


L 


004 16 


439 


490 C 


005 16 


233 




005 16 


121 


519 


004 16 


513 





005 16 


549 



855 0.162 



620 





782 


(J 


630 





080 


c 


84U 


(J 


530 


u 


062 


(J 




II 


927 






472 


0.270 


16 


374 


0.404 


SDSS 


760 


Null 




190 


Null 


CTIO 


792 


0.126 


15 


555 


Null 


SDSS 


850 


Null 


16 


750 


Null 


SDSS 


980 


0.230 




890 


Null 


SDSS 


611 


0.029 


n 


488 


0.033 


APO 




0.145 








APO 


071 


0.184 


16 


200 


Null 


APO 


009 


0.070 




876 




APO 


716 


0.152 


15 


530 


Null 


APO 



:opically Confirmed Sut 



SDSS J080823.63+171253.7 
SDSS J085039. 35+274536. 5 
SDSS J091018.78+292147.9 
SDSS J09 1033. 79+02021 1.2 
SDSS J 09 2 744. 74 +232402. 8 
LSPM J0943+6153 
SDSS J 100839. 65+064124. 2 
SDSS JlOlOll.01 + 175533.6 
SDSS J 1022 19. 30+052341. 6 
SDSS J 105007.88-002838.0 
SDSS J105100. 24+034120. 4 
SDSS J 113404.24+051234.4 
SDSS J 113458. 27+025400. 7 
SDSS J 115433. 70+094048. 7 
SDSS J 120405. 26+042145. 3 
SDSS J122216. 21+222739.0 
SDSS J 124425. 94-014425. 2 
SDSS J124930.69+030600.1 
SDSS J131014. 60+115508. 5 
LSPM J1412+0439 



08 50 39. 

09 10 18. 
09 10 33. 
09 27 44. 
09 43 38. 

08 39. 
10 11. 
22 19. 
50 07. 
51 00. 

34 04. 

34 58. 

54 33. 
2 04 05. 
2 22 16. 
2 44 25. 

2 49 30. 

3 10 14. 

2 29. 



62 +17 12 54. 

38 +27 45 36. 

81 +29 21 48. 

79 +02 02 11. 

76 +23 24 03. 

19 +61 53 22. 

64 +06 41 24. 

03 +17 55 34. 

30 +05 23 41. 

88 +00 28 38. 

25 +03 41 20. 

25 +05 12 34. 

28 +02 54 00. 

71 +09 40 48. 

27 +04 21 45. 

21 +22 27 39. 

95 -01 44 25. 

69 +03 06 00. 

60 +11 55 08. 

.00 +04 38 39. 



15 


760 





005 


4 


874 


15 


739 


t 


005 


4 


735 






t 




6 


069 


15 


859 


t 


005 


4 


902 


16 


786 


C 


008 


5 


880 






t 




5 




17 


262 


t 


009 


6 


267 




783 


t 




6 




16 


930 


t 


007 


5 


965 






t 




6 




16 


176 


C 


006 


5 


262 


17 


281 


t 


009 


6 


286 






t 


005 




595 


16 


542 


t 


006 


5 


424 




250 


t 






925 


IS 


245 


C 


014 


6 


965 


le 


503 


t 


017 


6 


900 




859 


t 




6 




n 


692 


t 


031 


7 


514 


16 


848 





007 


5 


902 



585 





003 14 


486 





004 14 


477 


463 


t 


004 14 


374 


i 


004 14 


385 




t 


004 15 


357 




004 15 


625 a 


524 


t 


003 14 


392 


t 


004 14 


326 


497 


C 


004 15 


378 


C 


004 15 


299 




t 


004 14 






004 14 


750 


881 


t 


004 15 


727 


t 


004 15 


382 




t 








005 16 


338 


560 


t 


003 15 


399 


t 


004 15 


347 


370 


t 


004 16 






004 16 


152 


939 


C 


004 15 


329 


c 


007 14 


767 


910 


t 


004 15 


751 


t 


004 15 


386 


055 


t 








008 13 


716 


993 


t 


004 14 


841 


t 


004 14 


781 




t 


004 14 






001 14 


251 


390 


C 


004 16 


176 


c 


004 16 


025 


320 


t 


004 16 


125 


t 


005 16 


D41 C 


809 


t 




585 








474 


t 


004 16 


J24 


i 


004 15 


777 


462 





003 15 


259 





003 15 


188 



c;tio 

CTIO 
CTIO 
CTIO 
CTIO 
APO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
CTIO 
APO 



CO 



. - SSS (Hambly et al. 2001), SDS (Mum 



I. 2004), LEP (Lepin 
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Deri\ 












WD 


Photometric 




Dist. 


Spec. 




Spectroscopic 




Name 


TefT (K) 


Comp. 


(pc) 


Type 


Teff (K) 


logs 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 




(7) 


(8) 


0312-084.. 


9080 ± 50 


He(+C) 


68.9 ± 11.2 


DQ 










c 


0351-002.. 


. . . 29,480 ± 260 


H 


129.5 ± 26.3 


DAI. 5 


36,960 ± 


190 


7.63 ± 0.03 


a 


0352-008 . . 


6160 ± 30 


H 


40.1 ± 6.7 


DA7.0 










b 


0412+065 . . 


.... 11,920 + 180 


H 


53.6 ± 9.5 


DA4.0 


13,720 ± 


280 


7.99 ± 0.03 




0811+250.. 


8220 ± 40 


He(+C) 


60.9 ± 9.9 


DQ 










d 


0840+243 . . 


5950 ± 30 


H 


38.0 ± 6.2 


DA7.5 










b 


0851+280.. 


7030 ± 40 


H 


63.5 ± 10.7 


DA6.0 


6780 


± 


90 


7.93 ± 0.17 




0900+204 . . 


7360 ± 40 


H 


52.2 ± 8.9 


DA6.0 


7190 


± 


40 


7.92 ± 0.06 




0907+221 . . 


5830 ± 80 


H 


40.3 ± 6.6 


DA7.5 










b 


0909+200 . . 


4950 ± 70 


H 


32.0 ± 5.3 


DA9.0 










b 


0920+012 . . 


6290 ± 60 


H(+Ca) 


33.4 ± 5.6 


DAZ 










e 


0921+315.. 


4810 ± 60 


He 


38.1 ± 7.3 


DC 












0948+192 . . 


5850 ± 100 


He 


51.5 ± 8.5 


DC 












1023+149 . . 


9190 + 90 


H 


83.1 ± 14.2 


DA5.5 


9180 


± 


50 


8.06 ± 0.06 




1025+197.. 


7330 ± 40 


H 


59.5 ± 10.2 


DA6.0 


7140 


± 


70 


8.34 ± 0.10 




1032+230.. 


6770 ± 20 


He(+Ca) 


42.2 ± 6.8 


DZ 










f 


1053+238 . . 


5110 + 110 


H 


39.4 ± 6.4 


DA9.0 










b 


1104+150.. 


6720 ± 110 


H 


37.8 ± 6.2 


DA6.5 


6710 


± 


90 


7.94 ± 0.20 




1110+292.. 


4810 ± 80 


He 


34.7 ± 7.1 


DC 












1116+288.. 


8850 ± 70 


H 


59.4 ± 10.2 


DA5.5 


8550 


± 


30 


8.37 ± 0.04 




1116-103.. 


7350 ± 40 


H 


42.8 ± 7.3 


DA6.0 


7360 


± 


50 


7.77 ± 0.09 




1119+289.. 


5690 ± 40 


H 


50.1 ± 8.2 


DA8.0 










b 


1134+209.. 


7860 ± 60 


H 


71.1 ± 12.3 


DA5.5 


7740 


± 


40 


8.01 ± 0.06 




1143+055.. 


6740 ± 40 


H 


43.1 ± 7.0 


DA6.5 


6840 


± 


60 


7.91 ± 0.10 




1145+304.. 


7060 ± 40 


He 


69.4 ± 11.5 


DC 












1151+246.. 


9060 ± 60 


H 


39.7 ± 6.6 


DA5.5 


8700 


± 


30 


8.74 ± 0.03 




1208+076 . . 


5430 ± 40 


H 


25.5 ± 4.4 


DA8.5 










b 


1224+321 . . 


6730 ± 40 


H 


38.7 ± 6.3 


DA6.5 


6530 


± 


80 


7.99 ± 0.15 




1229+151.. 


5840 ± 40 


H 


57.3 ± 9.4 


DA7.5 










b 


1239-072.. 


9970 ± 100 


H 


82.9 ± 13.7 


DA5.0 


10,100 ± 


50 


8.23 ± 0.05 




1257+185.. 


9500 ± 70 


H 


65.4 ± 11.1 


DA5.5 


9630 


± 


30 


8.50 ± 0.03 




1259+262 . . 


7050 ± 150 


H 


65.1 ± 10.9 


DA6.0 










b 


1307+143 . . 


9000 ± 40 


H 


55.8 ± 9.6 


DA5.5 


8610 


± 


30 


8.21 ± 0.04 




1322+092 . . 


8440 ± 40 


H 


56.3 ± 9.7 


DA5.5 


8210 


± 


40 


8.11 ± 0.07 




1338+052 . . 


4360 ± 50 


He 


13.7 + 2.7 


DC 












1401+069.. 


7140 ± 50 


He 


60.0 ± 9.9 


DC 












1402+065 . . 


... 26,190 + 210 


H 


116.6 ± 23.1 


DA2.0 


27,490 ± 


170 


7.83 ± 0.03 


a 


1404+163.. 


6810 ± 50 


H 


56.9 ± 9.2 


DA6.5 


6560 


± 


100 


7.65 ± 0.02 




1408+029 . . 


5570 ± 90 


H{+Ca) 


26.5 ± 4.3 


DAZ 










g 


1409+223 . . 


7160 ± 50 


He 


64.5 ± 10.7 


DC 












1419+062.. 


... 11,210 + 100 


H 


78.8 ± 13.6 


DA4.5 


11,620 ± 


110 


8.42 ± 0.04 




1425+057.. 


7060 ± 30 


He 


48.5 ± 8.1 


DC 












1434+159.. 


. . . 15,690 ± 170 


H 


112.8 ± 20.9 


DA3.0 


17,940 ± 


170 


8.06 ± 0.03 




1457+249 . . 


. . . 12,820 ± 90 


He 


85.6 ± 14.8 


DBA 


13,510 ± 


100 


8.09 ± 0.08 


h 


1516+128.. 


9290 ± 80 


H 


86.6 ± 14.6 


DA5.5 


9290 


± 


40 


8.12 ± 0.05 




1519+141.. 


... 10,160 + 60 


He 


78.6 ± 13.0 


DC 












1531+024.. 


8510 ± 60 


H 


49.1 ± 8.5 


DA5.5 


8220 


± 


30 


8.61 ± 0.04 




1547+481 . . 


6780 ± 40 


He 


51.9 ± 8.6 


DC 












1618+208.. 


8490 ± 60 


H 


61.3 ± 10.7 


DA5.5 


9070 


± 


30 


8.16 ± 0.04 




1630+089 . . 


5740 ± 20 


H 


13.2 ± 2.2 


DA8.0 










b 


2051-051.. 


. . . 10,470 ± 60 


He(+Ca) 


48.1 ± 8.0 


DZ 










' 


2102+233.. 


... 11,080 + 110 


H 


64.5 ± 10.5 


DA4.5 


12,040 ± 0.08 


8.36 ± 0.02 




2119+040.. 


5150 + 50 


H 


22.1 ± 3.6 


DA9.0 










b 


2144+156.. 


7730 ± 50 


H 


47.0 ± 8.1 


DA5.5 


8340 


± 


40 


8.49 ± 0.05 




^ T^ff from the 


SED fit is unreliable for th( 


hottest WDs 
















A combination 


of weak Balmcr lines and/ 


or noisy speet 


ra prohibited a ic 


iablc spcct 


ral fit. 










"^ Best fit model 


includes log (C/Ho) = -4.2 


± 0.2. 
















Best fit model 


includes log (C/He) = -5.1 


± 0.2. 
















** Best fit model 


includes log (Ca/H) = -9.0 


± 0.2. 
















f Best fit model 


includes log (Ca/Ho) = -9. 


) ± 0.2. 
















^ Best fit model 


includes log (Ca/H) = -9.0 


± 0.2. 
















^ Best fit model 


includes log (H/He) = -5.9 


± 0.1 
















' Best fit model 


includes log (Ca/He) = -10 


5 ± 0.2. 
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TABLE 3 
















Survey Summary 










Samp. 


PM 


Sky'^ 


All Objs. Selected'i # Obs. # Obs. 


# SDSS 


# Publ. 


# Publ. 


Remaining 




Source 


% 


N N WDs contain. 


Spcct. 


WDs 


OTHER 


N 


(1) 


(2) 


(3) 


(4) (5) (6) (7) 


(8) 


(9) 


(10) 


(11) 


Pass I'' 


USNO 


23 


97137 1028 19 18 


12 


105 


15 


859 


Pass 2^ 


LSPM 


29 


19408 190 6 2 


16 


71 


9 


86 


'^ Initial Coi 


istraints: PM 


> 0.10 a 


rcsec yr'l, Docl. < +30 deg.. g < 19.5, 2MASS mat, 


ch, SDSS DRS 


footprint 






*> Initial Co. 


nstraints: PM 


> 0.15 a 


rcscc yr^l, Dccl. > deg., g < 19.5, 2MASS match. 


SDSS DRS footprint 






^ Estimated 


by applying i 


nitial doc 


lination constraints to SDSS DR8 footprint. 










'^ Target sol 


ection criteria 


(see sect 


ion[2); RPM cut, J - Kg cut, model-adherence cut. 











Multi- Survey New White Dwarfs 



11 




4000 4500 5000 5500 6000 6500 



Wavelength (A) 



Ht Hi! H7 
I I 



H/3 
I 



tlO 



xi 

0) 



o 

2: 



j«lil/WMW»*W»****^r**'**»~--w*~>-»»_, 1104+150 

,1 , , 0907 + 221 




,««**---*'-*^ 





#«l##i*W^t**w^^ 



4000 4500 5000 5500 6000 6500 



Wavelength (A) 



Fig. 1. — Hydrogen-line DA WDs, plotted in order of decreasing T^g from top t o bot tom. Balmer lines are indicated and WD designations 
are labeled above the spectra. WD 1419+062 is included as discussed in Section 14.51 Spectra obtained at the CTIO 4m Blanco Telescope. 
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Fig. 2.— 
Telescope. 



DQ {top panel) and DC {bottom panel) WDs, plotted in order of decreasing T^g . Spectra obtained at the CTIO 4m Blanco 
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0920+012 




log (Ca/H) = -8.98 ± 0.20 
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Fig. 3. — DAZ WDs with Ca II H & K features. Model fits are displayed in the inset panels. Spectra obtained at the CTIO 4ni Blanco 
Telescope. 
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Fig. 4. — Contaminant spectra plotted in order of decreasing Tcfi . Spectra obtained at the CTIO 4in Blanco Telescope. 
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Fig. 5. — Spectra obtained from the ARC 3.5m Telescope at APO. Left panels represent the blue channel and right panels represent a 
subsection of the red channel centered on Ha. Top panels are new WD discoveries while the bottom panels are two contaminants, likely 
SDs. All spectra are plotted in order of decreasing T^ff . 
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Fig. 6. — Plot of select new DA WDs with spectroscopic determinations of log g and T^ff {filled stars) - from left to right: 0412+065, 
2102+233, 1419+062, 1239+072, 1257+185. Also plotted are WDs that are not observe d to vary (filled circles) and ZZ Ceti confirmations 
{open circles). The dashed lines represent the instability strip boundaries as defined bv lGianninas et al.l IJ2011I ). 
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Fig. 7. — Light curve (top panel) and FT (bottom panel) for WD2102+233. Object was observed from the CTIO 0.9m telescope at a 
cadence of ~ 50 seconds. From the FT, the dominant mode wras found to be ^ 800 seconds 
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Fig. 8.— 
cadence of 



2 2.5 3 

Frequency (mHz) 

Light curve (top panel) and FT (bottom panel) for WD1419+062. Object was observed from the CTIO 1.0m telescope at a 
~ 100 seconds. From the FT, no dominant mode was observed above the noise level out to the Nyquist frequency. 
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10 



'Hg= 15.136 + 2.727(ff - i) 




2.5 



Fig. 9. — Reduced proper motion diagram showing WDs from this study {open squares) and contaminants {filled squares). The dashed 
Hne represents the cut adopted from lKilic et aLl l|200g ) above which targets were discarded. Contaminants above the curve are discussed in 
Section [S] Modeled cooling curves for pure- hydrogen WDs with log g = 8 and Vtan = 30, 40, and 150 km s~^ are plotted as solid curves. 
Background stars were taken from SDSS for reference. 



